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Abstract 



We consider cosmological models of Bianchi type. In particular, we are interested 
in the a-limit dynamics near the Kasner circle of equilibria for Bianchi classes VIII 
and IX. They correspond to cosmological models close to the big-bang singularity. 

We prove the existence of a codimension-one family of solutions that limit, for 
t — )• — oo, onto a heteroclinic 3-cycle to the Kasner circle of equilibria. The theory 
extends to arbitrary heteroclinic chains that are uniformly bounded away from the 
three critical Taub points on the Kasner circle, in particular to all closed heteroclinic 
cycles of the Kasner map. 
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1 Introduction 

The Einstein field equations 

Ral3 — al3 = Tap 

restricted to spatially homogeneous, non-isotropic space-times gap with an ideal non- 
tilted fluid yield a system of ordinary differential equations, the Bianchi class-A model 
(EH) |WH89llWE05] . 

The a- limit, t — )■ — oo, of this system corresponds to the dynamics near the big-bang 
singularity. The dynamics in this limit, however, is not yet fully understood. 

It has been conjectured [Mis69t IBdSR86] that the dynamics follows the (formal) 
Kasner map deflned on the Kasner circle of equilibria and given by heteroclinic connections 
to the equilibria on the Kasner circle. Equilibria on the Kasner circle represent self- 
similarly expanding space-times. A trajectory close to a formal heteroclinic sequence 
thus corresponds to a space-time that is close to different self-similar spacetimes as it 
approaches the singularity in backward time-direction: a tumbling universe. 

At least for Bianchi class-IX solutions the Bianchi attractor formed by the union 
of the Kasner circle and its heteroclinic orbits has been proven to indeed be a (global) 
attractor for trajectories to generic initial data under the time-reversed flow [RinOl] . 

Rigorous results on the correspondence of the dynamics close to the attractor and the 
formal sequences of heteroclinic orbits to the Kasner circle are still missing. They would 
facilitate the discussion of the Belinskii-Khalatnikov-Lifshitz (BKL) conjecture of spatial 
decoupling close to singularities as well as the Misner hypothesis on the development of 
(spatially) homogeneous/isotropic space-times due to mixing near the initial singularity. 
See also |HU09] for a recent survey. 

In this paper we make the flrst step towards a rigorous description of the a-limit 
dynamics of the Bianchi system. We describe the set of initial conditions near the Bianchi 
attractor that follow the (up to equivariance) unique period-3 heteroclinic cycle. In fact 
we prove that this set forms a co dimension-one Lipschitz manifold, see theorem 14.21 

We note that this result does not depend on monotonicity arguments and therefore 
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covers both class- VIII and class-IX solutions approaching a period-3 heteroclinic cycle. 

We will start by reviewing the Bianchi system in section |2J in section 12] we study 
the passage near a line of equilibria in a generalized context. This yields a local map 
between sections to the period-3 heteroclinic cycle near the Kasner circle. In section H] 
this local map is combined with the global excursion given by the heteroclinic cycle. We 
obtain a return map with a fixed point representing the heteroclinic cycle. The stable 
manifold of this fixed point, i.e. the set of all solutions converging to the fixed point under 
iterations of the return map, represents all solutions following the heteroclinic cycle in the 
Bianchi system. We construct this (local) stable manifold as the limit object of a graph 
transformation. This yields the claimed Lipschitz set of trajectories with a-limit dynamics 



following the period-3 heteroclinic cycle, theorem 14.21 Finally, we discuss generalizations 
to this result as well as open problems in section |5l 

Generalizations include heteroclinic cycles of arbitrary period as well as non-periodic 
heteroclinic sequences that do not approach the singular Taub points. Matter models 
between dust and radiation can be included. 

Shortly after the submission of this article two more results on the same problem have 



appeared. While Beguin BeglO shows existence of solutions to non-periodic trajectories 
of the Kasner map that remain bounded away from the Taub points like in our treatment, 
Reiterer & Trubowitz jRTlO] construct solutions of the full system near trajectories of 
the Kasner map that pass arbitrarily close to the Taub points. However, their analysis 
seems to be restricted to the vacuum case. 



2 The Bianchi model 



We consider the Bianchi class-A model on (A^i, N2, N^, E+, E_) G with time-derivative 
' = d/dr. In terms of the spatial curvature variables Ni and the shear variables S± it 
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reads 



N[ = (g-4S+)iVi, 

= (g + 2S+ + 2v^S„)Ar2, 

iV^ = (g + 2S+-2v^S_)iV3, (2.1) 

= -(2-g)S+-35+, 

E'_ = -(2 -35_. 



The abbreviations 



q = 2(S2 +S2) + i(37-2)n 



2 ' 

= 1 - - S2 - ir. 



K = I (A^f + iV| + iV| - 2 (A^iA^s + iV2iV3 + iVsiVi)) , (2-2) 
5+ = 1 ((iV2 - iVa)' - iVi (2iVi - iV2 - iVa)) , 
5_ = '^VS{N3-N2){N,-N2-N^). 

include the deceleration parameter q, the density parameter Q, and the curvature param- 
eter K. The fixed parameter | < 7 < 2, given by the equation of state of an ideal fiuid, 
describes the uniformly distributed matter. For example, a value 7 = 1 represents dust, 
whereas 7 = 4/3 represents radiation. 

For a derivation of these equations from the Einstein field equations see [WH89j . 
|WE05j . the appendix to |Rin01] . or |Rin09] . |] 

The resulting fiow in Q yields 

n' = (2g- (37-2))fi. (2.3) 
The invariant set {Q = 0} corresponds to the 4-dimensional vacuum model 

N[ = 2(1 -ir-2E+)iVi, 

= 2(1 - + S+ + V3E„)iV2, 

iV^ = 2(1 - + S+ - V3S_)iV3, (2-4) 

S' = -2A'S, -35+, 



S'_ = -2A'S_-35_, resp. S_ = ± Jl - AT - S 



iNotc that |WE05| uses a slightly different scaling from |WH89j and }Rin01| . nI^^^ ^ ml^"^ with 
all other variables being the same. We use the scaling of |WH89) . 
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Bianchi Class 


Ni 


N2 




I 











II 


+ 








VIo 





+ 




VIIo 





+ 


+ 


VIII 




+ 


+ 


IX 


+ 


+ 


+ 



Table 2.1: Bianchi classes given by the signs of the spatial curvature variables iVj. Re- 
maining cases are related by equivariance. 

Symmetries are given by permutations of {A'^i, N2, N^} together with appropriate lin- 
ear transformation of S+, S_ corresponding to a representation of S3 on R,^. Together with 
the reflection {Ni, N2, N^) {—Ni,—N2,—N3), the system yields a 6*3 x Z2 symmetry 
group. 

Note the classification of restrictions of the dynamical system to the various invariant 
regions, see table 12.11 

The Kasner circle K, = {Ni = N2 = N3 = 0, Q = 0}, Bianchi class I, consists 
of equilibria. The attached half ellipsoids Tik = {N^ 7^ 0, A'; = = 0, ^2 = 0}, 
{k,l,m} = {1,2,3}, Bianchi class II, consist of heteroclinic orbits to equilibria on the 
Kasner circle, see figure [2?n The projections of the trajectories of Bianchi class-II vacuum 
solutions onto the S±-plane yield straight lines through the point (S+,S_) = (2,0) in 
the cap {A'^i 7^ 0, N2 = = 0}. The projections of the other caps are given by the 
equivariance. 

Away from the singular points, T^, k = 1,2, 3, the Kasner circle /C is normally 
hyperbolic with 2-dimensional center-stable manifold given by the family of incoming 
heteroclinic orbits. 

The Kasner map $ : /C — t- /C is defined as follows: for each point g+ G /C \ {Ti, T2, T3} 
there exists a Bianchi class-II vacuum heteroclinic orbit q{t) converging to g+ as t — )■ 00. 
This orbit is unique up to reflection (A''i, N2, N3) {—Ni, —N2, —As). Its unique a-limit 
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Figure 2.1: Heteroclinic caps of vacuum Bianchi II solutions to the Kasner circle. 
g_ defines the image of g+ under the Kasner map 



Including the three fixed points, $(Tfc) := T^, this construction yields a continuous map, 
$ : /C — 7- /C. In fact $ is a non- uniformly expanding map and its image $(/C) is a double 
cover of /C. 

The main goal are rigorous results on the correspondence of iterations of the Kasner 
map $ to the dynamics of nearby trajectories to the Bianchi system (12.11) with reversed 
time, i.e. in the a-limit t — — oo. 

There exists a 3-cycle of heteroclinic orbits, i.e. a fixed point of unique up to 
equivariance. In the following we choose the cycle given by a heteroclinic orbit q{t) in the 
Bianchi class-II vacuum cap {A^^i > 0, = ^3 = 0} with S_ > and its images under 
the equivariances given by cyclic permutation of Ni and rotation by 27r/3 in (S_|_, see 
figure 12. 2[ The a-limit of q{t) is given by 



After factoring out the equivariance, the heteroclinic orbit q(t) becomes a homoclinic orbit 
on the orbit space of the equivariance. We will therefore discuss the dynamics of nearby 
trajectories by studying the return map to a transverse cross section to h{t). This is 
done in section |H 



$(g+) := g. 



(2.5) 




(2.6) 
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The linearization of (12.1 p at the Kasner circle yields 



2 - 4S+ 
2 + 2S+ 








2v^S_ 
2 + 2S+ 








2V3S_ 

3(2 -7)S^ 
3(2-7)S+S_ 








3(2-7)S+S_ 
3(2-7)S^_ 
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We find eigenvalues 

/xi = 2-4S+, 

/i2 = 2 + 2S+ + 2v^S_, 

/i3 = 2 + 2S+-2V3S_ 
to eigenvectors d^i, dN2, tangential to the Bianchi class-II vacuum heterochnics. 
Additionally, there is the trivial eigenvalue zero to the eigenvector — + T,^d^_ 
tangential to the Kasner circle /C. The fifth eigenvalue /i^ = 3(2 — 7) > corresponds to 
the eigenvector + transverse to the vacuum boundary {Q = 0}. 

At the Kasner equilibrium g„ of the 3-cycle (12.61) we have 
3 3 

Note that both unstable eigenvalues are stronger than the stable one, < —/is < fi2 < f^i, 
and that the heteroclinic orbit belonging to the 3-cycle is tangent to the strong unstable 
direction d^^^. 

In fact at every point on /C \ {Ti,T2,T3} there is one negative eigenvalue and it is 
weaker than the other two positive eigenvalues among fii, fi2, fi-s- At the singular Taub 
points {Ti,T2,T3} two eigenvalues change their signs simultaneously. 

We therefore study the local passage near such nonsingular points in general systems 
in the following section, before combining it with the global excursion map and applying 
the results to the particular Bianchi system in section |H 



3 Local map 

In this section we study the passage of trajectories under a general flow near a line of 
equilibria with eigenvalue constraints (13. 4p consistent with the Kasner circle in the Bianchi 
system. We will collect estimates on expansion and contraction rates to establish Lipschitz 
properties of the local map between sections to a reference orbit given by the passage near 
the line of equilibria, see theorem 13.91 at the end of this section. 

Consider a vector fleld, k > 4, 

x' = fix), X e R^ (3.1) 
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Figure 3.1: Local passage : S'"" ^ S 



■'out 



near a point (w.l.o.g. the origin) on an equilibrium line 



/ (0,0,0, =0, Vx.eR. 



(3.2) 



Assume that the linearization at the origin (and then by continuity locally all along the 
line) has the form 

^ /in(Xc) 

-fIs{Xc) 

-/iss(Xc 

* * * 



\ 









(3.3) 



with 



(3.4) 



< /i„ < /is < fiss- 

We denote x = {xu,Xs,Xss,Xc)^ according to the above splitting. We also abbreviate the 
vector of stable components as Xs,ss = {xsiXss)"^- 

The aim is to study a local map from an in-section = {xss = s} to an out-section 
^out _ _ ^ j ^ figuj-e 13.11 This corresponds to the passage near the 

Kasner circle in the Bianchi system in backwards time direction. (We reversed the time 
direction to obtain a well defined local map.) 

Rescaling coordinates yields a system 



x' = A{xc)x + eg{x), 



(3.5) 
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with e arbitrarily fixed and g quadratic in {xu,Xs,Xss)- The local map 

is given by the first intersection of the solution of fl3.5p to the initial value = 
with the out-section {x„ = 1}. The singular points in the intersection of the 
stable manifold of the equilibrium line with the in-section are mapped to the respective 
points in the intersection of the unstable manifold of the equilibrium line with the out- 
section. Time at the in-section is usually set to zero. Time at the out-section depends on 
the initial condition and is denoted by 

For any Xc fixed and any /c G IN, the equilibrium (0, 0, 0, Xc)""" possesses a one- 
dimensional unstable manifold W'^{xc) and a stable manifold iy^'*'^(xc). 

Since the stable and unstable manifolds as well as the strong stable foliation of the 
stable manifold are they can be fiattened, see e.g. |SSTC98] . Theorem 5.8. By a 
change of coordinate the stable / strong stable / unstable manifolds to the equilibria 
locally coincide with the respective eigenspaces, in particular the following subspaces 
become invariant: 

W{xc) = {xs = Xss = 0,Xc^xed}, 
W''''{x^) = {x„ = 0,Xe fixed}, (3.7) 
W'^lxc) = {xu = Xs = 0, Xc fixed}. 

In fact, in the Bianchi system, W^{xc) and W^^{xc) coincide with the class-II caps formed 
by families of incoming and outgoing heteroclinic orbits. 

The local map is well-defined on the in-section 

S'" = { (x|r, xT, xZ, X-) I X- = 1, < x|r < 1, |xr I < 1, K\ < 1 } (3.8) 

see lemma ESI below. The singular points of the local map thus become the set {xj^ = 0} 
and we define: 

vl/i°=(xL" = 0,x^xn = (xr,x:f,xr) := (0,0,x-). (3.9) 
Let Pu, Ps, Pss, Ps,ss '■= Ps + Pss, Pc be the eigenprojections with respect to A{xc)- 
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Then (13. 3p yields 



A{xc) = Df ({0,0,0, X, 



' IJ^uiXc) 






^ 

fls{Xc) 

-fissile) 





\ 
and due to fl3.7l) the higher order terms have the form 

Pug{x) = XuCjuix) e R 

gsi{,x) gs2{,x) 



(3.10) 



~ Xu9u\X) 
Ps,ssg{x^ Q s,ss{x)x s^ss 

Pcdi.x) Xu9cix^X ggg 



Xii 



gssi[x) gss2[x) 



gci[x) gc2[X) 



e R, 



with C -functions gu, gs.ss, vanishing along the line of equilibria, and C -function g 



(3.12) 



In particular 

|^ci(a;)|, |^c2(a;)| < C, 
\gu{x)\, \gsi{x)\, \gs2{x)\, \gssi{x)\, \gss2{x)\ < Cmax(|x„|, \xs\, \xgs\) 

for some constant C > independent of e and x & where U is some local neighborhood 
of the origin. Similarly, g satisfies Lipschitz bounds 

\gc{x) - gc{x)\, \gu{x) - gu{x)\, \\gs,ss{x) - gs,ss{x)\\ < C\\x - x\\. (3.13) 

Norms of vectors are always taken as £i-norms, e.g. 



\x\ 



\Xc\ ~t~ \Xu\ ~t~ \\Xs 



\Xc\ ~t~ \Xu\ ~t~ \Xs\ ~t~ I'^ssr 



We choose 



U = {-2,2f. (3.14) 
All further estimates will use this rescaled system (13.51) with flattened invariant manifolds 



(13. 7p in the local neighborhood U. They will be valid for all e < Eq and suitably chosen 
Eq. In the original system (13. ip Eq bounds the size of the neighborhood of the origin in 
which this local analysis is valid. 

Remark 3.1 The invariance o/iy**(xc) implies that gs2{xu = 0,Xs,Xss,Xc) = 0. We do 
however not exploit this fact in our analysis below. 
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Proposition 3.2 Let 



be the eigenvalues of l{3.10\) at the origin. Then for all < a < 1 there exists an Eq > 
such that for all e < Eq in 113. 5\} and x ElA 



fluiXc) fis{Xc) fiss{Xc) , __i 

a < , , < a 

Proof. The distinct eigenvalues in the linearization of the original system f l3.3p depend 
different iably on Xc, as long as f l3.4p holds. For the rescaled system f l3.5p with small Eq 
this provides bounds in U: Indeed, there exists a constant C > independent of Eq, e, 
such that 



d 
dxc 



fig 



< eC. 



(3.15) 

CXI 



The scalar function 6{x) := /i^ (//^(xc) + Egu{x)) is therefore C'^ ^ and close to 1. 
The vector field 



x' = 9{x)f{x) 



-fix) 



has the same trajectories as the original vector field and all previous considerations remain 
valid. Thus we can assume, without loss of generality, that 6{x) = 1 inU, i.e. 



l^uyXc 



gu{x) = 0. 



(3.16) 



At this step we have made use of the fact that the origin possesses exactly one unstable 
eigenvalue. The vector field to consider then has the form 



X 



Xu 

I 

s,ss 

x' 



f^u Xu^ 

^'S,ss{Xc) Xs^ss ~l~ ^ gs,ss(^x) Xs^ssy 
^ Xu gdx) Xg^gg. 



(3.17) 



Note again the abbreviations 



X Q 



X s 



^ flg{Xc) 
^ f^ggiXc) 



gs,ss{x) 



gsi{x) gg2{x) 
ggsi{x) ggg2{x) 



9c[x) 



gci[x) gc2[x) , 
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to facilitate later generalizations to higher dimensions of the stable component Xs,ss- In 
particular note the diagonal form of the linear part ^s,ss- 

Lemma 3.3 Define 

Then for all < a < 1 there exists an Eq > such that for all s < Sq and x & U 



unt 



der the dynamics 1^3.11 ). In particular, with x(0) = on the in-section as long 



as the trajectory remains inside U 

X{t) < exp(-a(^, - ||x"^J| < 2exp(-a(/i, -^„)t) 

Remark 3.4 The quantity x describes the "distance" from the critical heteroclinic cycle. 
More precisely, x is related to the distance from the set IJx-c {W''''^'^{xc) U W^^Xc))- Close 
to the critical heteroclinic cycle this set is however under the return map squeezed onto the 
critical cycle, so for large times x measures in fact a "distance" to the critical heteroclinic 
cycle. 

Proof. Vectorfield (13.171) yields 

X < -Mxc) - |J,u)x + '^£\\9s,ssix)\\x■ 
The function g\ii is bounded uniformly in e, see fl3.12p . Additionally, proposition 13.21 
provides bounds —fis{xc) < —ofig for arbitrary < <5 < 1 and Eq small enough. By a 
suitable choice of a between a and 1 we estimate 

X < -a{^s 

Integrating this differential inequality and using the initial value yield the second claim. 

CXI 

We need estimates for the passage time as well as bounds on \xc\, \xu\, and 



These are obtained in the following lemmata and summarized in corollary 13.71 
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Lemma 3.5 For all < a < 1 there exists an Eq > such that for all e < Eq and 
x{0) = X™ in the in-section S™, see ^3. <S]) . as long as x{t) remains in U under the flow to 



the vector field \3.1'T\j 



x^{t) = exp{fij)x'^, (3.18) 

= \xs{t)\ + \xssit)\ < exp(-a/i,t) ||x!,'^^J|, (3.19) 

2fC 

K{t)-x^\ < — (3.20) 

Here, C is the uniform (in x and e) hound on g\u from 113. 

Proof. The unstable component (13.181) is given directly by the vectorfield. 
On the stable component, we obtain the estimate 

■f^{\Xs\ + \Xss\) < {-fis{Xc) +2e\\gs,ss{x)\\) {\Xs\ + \Xss\) ■ 

The function g\K is bounded uniformly in e, see (I3.12p . Together with the estimate on 
the eigenvalues of proposition 13. 2[ this yields the bounds (I3.19P on the stable component. 

The center component of (I3.17P is estimated by 



f P,g{x{s)) ds < e f'\P,g{x{s))\ds 

2 



I c I — ^ 

< eC I \x{s)\ds < eC—^ 

Here we used the estimate of lemma 13.31 on x- ^ 

Corollary 3.6 There exists an Eq > such that for all e < Eq and x{0) = in the in- 
section see Ii3.8\) . the trajectory x{t) under the flow to the vector field Iji3.n\ ) remains 



inU as long as \xu\ < 1, i-e. all along the passage defining the local map The passage 
time is given by 



f^u I I 

Proof. We choose Eq smaller than ^a{fis—f^u), see lemma 1331 Then trajectories starting 
in S'"^ cannot leave U unless x„ becomes larger than 1, see (13.201) . (13.191) . Furthermore, 
(I3.18P ensures that Xu must grow beyond 1. Thus every trajectory starting in intersects 
the out-section {xu = 1} before leaving U. Setting Xu(t^°^) = 1 in (I3.18P determines the 
passage time t^"^. [xi 
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Corollary 3.7 The local map l[3.6[\3.[^} . i.e. the local passage on the closed in-section 
including the singular line {x™ = 0}, is continuous. For a// < a < 1 there exists an 
> such that for all e < Eq the following estimates hold 



ixr^-xi"! < — -\x. 



u I ' 



a (/is - 

||^out|| _ |™Out| I I outi ^ \ \n\aiis/nu-l\^m\\\ in ii ^ |^in|Q!Ats/Mi/-l/'|^in| , I m|\ 

l-^S I "T l-^SS I — \-^U I l-^M I — \-^U I \\-^ U I S \)- 

Thus the drift along the line of equilibria is arbitrarily small and the distance from the 
orbit to the union of the stable and unstable manifolds shrinks arbitrarily fast, close to the 
critical orbit. 

Proof. The estimates follow directly from lemma 13.51 applied to the local passage time 
given by corollary 13.61 They also establish continuity of the local map at the singular 
line {x™ = 0}. Note < x™ < 1 and = 1 on the in-section. Note further that 
P = afis/ ftu — 1 > for a chosen close enough to 1. ixi 

To obtain Lipschitz-bounds for the local map \lf^°'^, we have to improve the estimates 
considerably. Consider two trajectories, x{t) and x{t), starting in the in-section S™. If 
x^ = then x°"*^ = 0, x°^* = x™, and the Lipschitz estimates follow from corollary 13.71 
Choose, without loss of generality, 

< x|r < X-. 

In order to compare the two trajectories, we want Xm(0) = x" = x„(0) and therefore start 
the trajectory x at negative time 

x(-tO) = x'", P := — In^. (3.21) 

Then Xu{t) = Xu{t) for all t > 0. In particular, both trajectories meet the out-section at 
the same time t^°'^ = — In v-j^ , see corollary 13.61 The next lemma gives some estimates 
for the initial segment of the trajectory x up to t = 0. 

Lemma 3.8 Consider ( [5*. i7p with initial values x(0) = x™ and x{—t^) = x™ in S™. 
Assume without loss of generality < x™ < xj^. There exists an Eq > such that for all 
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e < Eq the following estimates hold: 



1, 



|x,(0)-xi"| + |x,,(0)-4"J < Ls 



I in _ -in I 



with suitably chosen constants Lg > 0, Lc > 0, uniform in x™, x'" and e. 



Proof. The Xc-component is estimated similar to lemma [3T51 

|xc(t)— x'"| < eC / |x(s)|ds < eC / exp(— — /i„)s)|x™| ||x 



Again C is the uniform bound on qIk. The last inequality uses the estimate of lemma 1X51 
Plugging in (13.211) and < 2 we obtain 



\x,{0)-x[^\ < 



2eC 



< 



< 



2eC 



(l - exp{-a{i2s - At«)t°)) 



2eCL 



1 - ^ . 



The last inequality with L := max{(/is — fiu)/f^u, 1} uses the general Lipschitz estimate 
1 — < maxj??, 1}(1 — ^) for arbitrary < < 1 and > 0. Finally, the choice 
Lc := 2CLa~^{fis — f'u)^^ yields the claim on the Xc-component. 

The Xs.ss-component is contracted at most with the rate of the strong stable eigen- 
value figs, perturbed by small higher order terms. Indeed, due to the diagonal form of the 
linear part, the variation-of-constant formula reads 



Xs,ss (0) 



V 



exp/°^-o -/is(xc(s)) ds 







\ 







exp/°j(, -Hss{xc{s)) ds 



X™ 
•^s.ss 



+ e 



V 



exp/° -fis{xc{cr)) da 











exp J° -fiss{xc{(T))da 



9s,ssi^x(^s^^Xg gg ds. 



We have already established uniform bounds on the nonlinearity, \g\ < C, and the con- 
traction rates, a fig < fis{xc) < fiss{xc) < -fiss, with < a < 1 arbitrarily close to 1. 
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Thus 



^(0)-x^%|| < (1 -exp(-i/i,,t°)j + eCy exp(-a/i,s) ds 

eC 

< 2 (l - exp(-i/i,iO)) + (l - exp(-a/ij°)) 



~in\ Mas/a^iiX / / ^in \ /Mu 




2 1 - 

(~in N 
1- ^ 

The last inequahty again uses the general estimate 1 — .^'^ < max{^!^, 1}(1 — for arbitrary 
< ^ < 1 and ^ > 0. The constants Li and L2 are independent of e and therefore yield 
the claim. ixi 

Theorem 3.9 The local map \E'^°'^ (EI3 \3.9\) . i.e. the local passage on including the 
singular line = 0, is Lipschitz continuous. Let 

||~in _ inll _ | -in _ inl , |~in _ inl , |~in _ inl 
IP II l-^M -^u I ~r -^c I ~r l-^s -^s I 

be the norm in the in-section S'"^. Then there exists an Eq > such that for all s < Eq 
and < < the following estimates hold: 

||™out ™out II |™out ™out I I |~out ™out I ^ T |™in|/9||~in t'^H 

11-^5,5^ -^SjSsll l-^s -^s I ~r l-^ss -^ss I — ^s\->^u I 11-^ In 

wi/i Lc > 0, L^, > 0, /3 > uniform in x, e. In particular, Ls\x^^\^ is arbitrarily small 
for a;™ small enough. 

The drift in the center direction can be made arbitrarily small by choosing a suffi- 
ciently small local neighborhood. The contraction in the transverse directions is arbitrarily 
strong by restricting the in-section to the part close to the primary object, i.e. the stable 
manifold of the origin. 

Proof. If X™ = then = 0, = x™, and the Lipschitz estimates follow from 
corollary 13. 71 Therefore let < x™ < x^^. Choose x(— = x™ according to fl3.2ip . Then 



x(0) is estimated by lemma l3^ and x(0) = x™. 

We abbreviate 

dc{t) := f (xe(s) - Xe(s))' ds, 4(0) = 



(3.22) 
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In particular, using lemma I3.8[ we estimate the initial value, 

4(0) < Lo^K-S^\ + ^S^,,-xZs^ < (3.23) 

as well as the distance in the center direction, 

\icit) - Xc{t)\ < dc{t) + \xc{fi)-x'^\ + \x';^ -x';^\ < dcit) + L2\\5^'' - ^^''W. (3.24) 

The constants Li, L2 > are uniform in x, x, e. 
Vector field f l3.17p yields the estimate 

(•^SjSs -^s^ss) f^s,ss(^'^c) (^■^ S,SS -^SjSs) ~l~ R 

with remainder term 

||-^|| — \\fJ'S,ssi^-^c) A''S,ss('^c) II ll-^SiSsll 

^ llfi'SjSS ('^) II II *^S,SS "^SjSS II ~l~ ^ II S'SjSS (2') fi'SjSS ('^) II II -^SiSS II 

— sC\Xq Xq\ II-'^s^ssII ~t~ ^^||^s,ss -^SjSsll £^C^||x x|| 11x^^5^11 

^^(|-fc -^cl ||^S,Ss|| ~1~ 11*^8,^8 *^S,Ss||)* 

Here, C is the uniform bound on g and uniform Lipschitz-bound on g and ^s,ss, see f l3.12p . 
fl37[3|) and f l3J[5|) . 

In the last step we have used 

rf' nr> n/-> rp I I rp /y» | 

II II I c '^cl ~r W'^s^ss ^s,ss\\ 

I = due to our choice fl3.21l) . 

For arbitrary < a < 1 we have a fig < /is(xc) < f^ssi^c) provided Eq is chosen 
accordingly. Hence 

dsity < \\Xs,ssit) - Xs,ssit)\\' 

< -afisds{t) + eC\xc{t) - Xc{t) \ \\xs,ss{t)\\. 
Analogous Lipschitz estimates on cjc yield 

deity = \xe{ty - xe{ty\ 

^ £ |Xji I II (3^) II II -^SjSS •^S,SS II ~l~ ^ l-^M I II (3^) QciS^^ II II -^SjSS II 

< eC\Xu{ty {\\Xs,ss{t) - Xs,ss{t)\\ + \Xc{t) - Xc{ty ||x,,3,(t)||) 

= eC expifijyx^l {ds{t) + \xc{t) - Xc(t)| ||a;,,s,(t) ||) . 
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Using estimate f l3.24p . the exponential bounds on Xs,ss of lemma [375| and < 2, we 

obtain 



ds{t)' < -afisdsit) + 2eC exp{-anst) [dc{t) + L2\\x''' - x' 
d,{ty < eCexp{fiJ)\x'^\ds{t) (3.25) 
+ 2eC exp{-a{ns - + LsHx^"^ - x' 



The Gronwall estimate of exp{afist)ds(t), 

ds{t) < exp{-afist) (^4(0) + 2eCtL2||x'" - + 2eC dc{s) ds^ , (3.26) 

can be plugged into the estimate of d^- 

dc{ty < 2eC exp{-a{fis - fJ^u)t)\x'^\ {dc{t) + eC d,{s) ds 

+ 4(0) + {eCt + 1)L2||£'" - 

Let — be an upper bound on dc, i.e. assume 

llx^^-x^^ll > sup d,{s). (3.27) 

0<s<t 

As dc{0) = 0, this assumption is valid on an initial segment. Then 

deity < 2£Cexp(-«(/i, - fiu)t)K\ {dcit) + 4(0) + (eCt + 1)L3||5'° - x'^) , 
with L3 = L2 + 1. Gronwall inequality yields 

dc{t) < ^{t) + ^{s)r]{s) exp (^J\{(r) dcr^ ds 

with 

rj{t) = 2eC exp{-a{fis - fJ'u)t)\x'^\ 

< 2eC exp{-a{fis - fJ'u)t) , 

pt 2£C 

/ r]{a)da < -, < s < t, 

e(t) = 2eC\x'^\ J' exp{-a{fis- fiu)s)[d,{0) + {eCs + l)L3\\x''' -x'^'W) ds, 

^ ^^H^KI (^^(0) + ^3||5^- - x-ll) + e'^^,^g^\x^:\L4x^- - x-ll 

< eL4|xjf| (4(0) + 

< eL4(Li + l)||i'° 
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The second last inequality needs Eq chosen small enough, L4 is a uniform constant. The 
last inequality uses the initial estimate fl3.23p and < 1. Finally 

dc{t) < m + exp (^^* r]{a) da^ ds 

< eL,{L, + l)\\x'^-x'^\\ ^l + 2eCexp{^^^^^)j\xp{-a{fi,-fi^)s)ds^ 

< eL^Wi:''' - x'^'W, 

(3.28) 

with L5 independent of e. Choose Eq small enough. Then eL^ <^ 1 and assumption fl3.27p 
remains valid up to the out-section. This yields the first claim of the theorem. 

Estimates f l3.28p and fl3.26p combine to 

< exp(-a/i,t) (4(0) + etLellx'" - 
Thus 

ci°"t _ 1 ) 

< \x':r^/^- ( L^tl^l^ + £:^|Uin _ ^inn 1^ 1 \ 
— \ 'rin // Tin / 

Note that /3 = afis/ yU« — 1 > for a chosen close enough to 1. This finishes the proof, ixi 



4 The return map 

In this section we define a return map for trajectories near the primary heteroclinic orbit 
q{t) considered as a homoclinic orbit relative to the equivariance. In this way a map on the 
in-section is defined which induces a transformation of Lipschitz-curves on this section. 
We first prove Lipschitz- and cone properties of the return map. The homoclinic orbit 
q{t) corresponds to a fixed point of this map. The stable set of the equilibrium under the 
return map is then constructed as the fixed point of the associated graph transformation. 
This yields the set of trajectories with a-limit dynamics following the period-3 heteroclinic 
cycle resp. the relative homoclinic orbit q{t). 

Locally near g„ we can use the general result on the local passage of the previous 
section and apply it to the Bianchi system. The time-reversed Bianchi flow satisfies 
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assumptions (13. 2[ 13. 3p . The out-section is then diffeomorphically mapped back to the 
in-section by the time-reversed (global) Bianchi flow after reduction to the orbit space of 
the equivariance group, 

^'^out out ovLt\ , ( in in m\ _ ^r^glob/ out out out\ (a i\ 

In particular the line {x°^^ = = 0} is mapped to the line {x™ = = 0} by the 
Kasner map. Thus expands the Xc-coordinate uniformly in a neighborhood of the 

3-periodic cycle q(t). 

Combining both maps (the local passage and the global excursion) we obtain a return 

map 

^ := ^Sloh Q ^loc . ^in ^ ^in^ ^^^^^ 

see also figure 14. 1[ Note, that due to expansion of the center component Xc under the 
return map, the domain of definition of this map is in fact a subset of S™. 

The heteroclinic cycle q(t) yields a homoclinic orbit after reduction to the orbit space 
of the equivariance group. It constitutes the unique fixed point (0, 0, 0) of the return 
map The Kasner circle and the attached heteroclinic cap is represented by the set 
{x™ = x^° = 0}. The return map becomes singular on the stable manifold of the Kasner 
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circle. In particular, we define on this singular set 

v^(0,a;r,xn = (0,0,$(4-)), (4.3) 

where $ is the Kasner map. 

Note that the Bianchi system (12.11) preserves the signs of A'^i, ^3, ^- Although the 
local analysis of the previous section did not use this special structure, it is compatible 
with it. To study Bianchi-IX solutions close to the 3-cycle of heteroclinics we therefore 
consider only the part of the in-section with x„ > 0, > 0. This part is mapped to the 
in-section, again with Xu > 0, > 0. 

In complete analogy one can treat Bianchi- VIII solutions close to a 3-cycle, although 
this 3-cycle is then composed of orbits inside heteroclinic caps of different sign. Indeed, 
each heteroclinic orbit to the Kasner circle has a mirror image under sign reversal of N^. 
Each of the 3 heteroclinic orbits of the cycle can be given a sign, the choice of signs 
determines whether the cycle lies in the boundary of Bianchi- VIII or Bianchi-IX domains. 

In fact, looking at a section to a given heteroclinic orbit of the cycle, thus fixing the 
sign of the first orbit, we have 4 quadrants in this section corresponding to the 4 remaining 
choices of signs of the 2 remaining orbits. In each of the quadrants, theorem 14.21 will yield 
a local Lipschitz manifold. 

Although the analysis of this section also holds for a 3-cycle in the boundary of 
Bianchi- VIII domains, we keep the notation corresponding to Bianchi-IX, i.e. all com- 
ponents Xs,Xss,Xu are and remain non-negative. These components correspond to the 
variables Ni,N2,N3 in the Bianchi system, the identification permutes along the cycle. 
The component Xc correspond to the variables T,±. 

In the following, we first prove suitable cone conditions for the return map \l/ and then 
apply almost standard graph transformation to obtain the claimed Lipschitz manifold. 
The nonstandard part of that graph transformation includes the singular line xj^ = of 
infinite contraction. 

Note the order of fixing the rescaling parameters: First Eq resp. e is fixed small 
enough to yield our estimates of the local passage with small Lipschitz constants, in 
particular eLc ^ 1. This amounts to a choice of the sections and E°^* in the original 
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(unsealed) eoordinates and also fixes the global excursion map 

^giob rpj^g^ ^ sufficiently 
small upper bound for is chosen, i.e. is restricted to a smaller section 

= {xeT''\0<Xu<6, 0<Xs<6, \x,\ < 6} (4.4) 

This makes the contraction of the local passage as strong as we like without changing 
v|/i°c^ vl/si°t). It also ensures that trajectories of interest stay close to the Kasner caps 
of heteroclinic orbits, thus we ensure that the global excursion vl>si°ti ^j^g domain of 
interest is as close to the Kasner map as we like. It also ensures that all non-singular 
trajectories in this domain indeed return to the in-section S™. 

Lemma 4.1 The return map ^4-2^ is Lipschitz continuous. Furthermore, there exist 
e > 0, 5 > 0, < 0" < 1, Ku,s > Ij (ind Kc > (1 — o"^)"^ > 1, such that the following cone 
conditions hold for * = o : S*" ^ S"\ 

Here is the closed in-section Ii3. 8\) corresponding to the choice of e, and a 
suitable subset lji4-4^ - 

The cones are defined for x G as 

~ {x ^ I \\Xu,s ~ Xu,s\\ ^ cr\xc ~ Xc\}- , / , ^\ 

^x'^ ~ {x ^ I \xc — Xc\ < cr||Xu^<j — x^i^sll}. 
The cone conditions are 

(i) Invanance: ^(C^) n t''' C (int C^J U {^x} and ^"^Cfe) H S'" C (int C^'') U {x}; 

(a) Contraction & Expansion: For all x G we have expansion in the center direction: 
|(\E'x)c — (^E'x)cl > Kc\xc — Xc\ and for all G Cj^'^ we have contraction in the 
transverse directions: \\xu,s — Xu,s\\ > Ku,s\\{'^x)u,s — {'^x)u,s\\- 

They hold for all, x, x, G S™. See also figure 

Proof. Lipschitz continuity of the return map follows directly from Lipschitz continuity 
of the local passage \Ef'°'^, see theorem 13. 9 [ as the global excursion v[/gi°t> jg gmooth. 

The cone conditions require the expansion in Xc-direction along the Kasner circle 
induced by the Kasner map. On the singular line, the global excursion is given by the 
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^u,s Expansion 
of Xc 
Contraction 

of Xqi_C 




Figure 4.2: Cone properties of the return map 



Kasner map $, 



Indeed, the local passage does not change x^ on the singular line, so this notation is 
consistent with f l4.3p . Therefore we can write 

with a smooth (3 x 2)-matrix \E'^^°*^(xs, x^s, Xc) and vector Xg.ss- We obtain the following 
Lipschitz estimates 



'^s,ss ~l~ ll-^ -^ll ll-^SiSsll) ) 



(x) (Xs^^s Xg^ss) 

lob 



< cs'°^riix 



with uniform bound (7^^°^ > HvE^s^"*^]!, ||Z}\['S'°^|| on the sup-norm of the higher order terms 
of the global excursion map and their derivatives. Note that x, x lie on the out-section 
and vl^s^"^ is a diffeomorphism from the out-section to its image on the in-section. 

Using \E'^°'^(x) and \&^°^(x) instead of x and x we get a similar estimate for the return 
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map ^(x) = 

< 0^'°'° {\\^'°'{x)s,ss - + - (4.6) 

< Csi°b - x|| + (eL, + 1 + L,|x„|'^)||x - x|| 

^ (^returnl _ ^11 

The second last inequality uses the estimates of the local passage of corollary 13.71 and 
theorem 13.91 for the choice (w.l.o.g.) < x^, < Xu- Note that the estimates of theorem 13.91 
are used in the form 

|*^°^(x)e-^'°'(x),| < eL,\\i-x\\ + \xc-x,\<{eL, + l)\\i-xl 

< Ls\xuf\\x-x\\. 

The constant (7'''^*™ is uniform in x, x in the in-section. Because (3 > 0, we have an 
arbitrarily strong contraction on S™, i.e. for x^ < 5, if we choose S small enough. 

The Kasner map $ in the original Bianchi system is expanding everywhere except at 
the Taub points Ti,T2,T3. In particular it is expanding on our in-section: 

|$(a)-$(6)| > Kc\a-bl 

for some uniform constant Kc > 1. 

Now choose with 1 < < and a with < cr < 1 such that Kc{l — a^) > 1. 
(The last relation is needed to obtain a contraction in theorem 14.21 ) 

Consider the cone in center direction with opening i} > 0, i.e. .,|| < '&\xc—Xc\- 

Then fl4.6l) using the local Lipschitz estimate of theorem 13.91 yield 

1(^(5) - ^(x))c| > |$(^'°^(x)c)-$(^^°'(x)c)|-C'°*™|x„|^||x-x|| 

> i^c|*'°=(x)c - *^°'(x)c| - C"-^*"™|x„|^||x - x|| 

> Kclxc — Xc\ — eLcKc\\x — x\\ — C'''^*"™|x«|^||x — a:|| 

> (J?e - {eL.K, + (1 + ^9)) \x, - x,\. 

For e and 6 chosen small enough, using Xu < S, we can achieve 



(4.7) 



< K,-[eL,K, + C'''''''^\xuf){l + l/a), 
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yielding the expansion not only in the cone C^, with -d = a < 1, but also outside the cone 
Cl'\ with § = 1/a. 

Furthermore, using again (14. 6p . we see the invariance of the cones. Indeed, assume 
again \\xu,s — Xu,s\\ < ''^\xc — Xc\, then we have 

\\'^(x) —"^(x) II < C"'''*"™^ l^llx — xll 

< C"^^t"'-=^|x«|^(l + i^)K-^\^{x)c - ^{x)c\. 

The last inequality uses the expansion in Xc, thus it is valid for -i? < 1/a. We choose 6 
small enough such that C"'°*"™|xti|^i^~^ < a/{l + a). Due to the monotone increase of 
^9/(1 + 1^) we also have C"'^'''''\xu\^K-^ < ^/{l + ^9) for all ^ > a. Thus we obtain the 
cone invariance 

||^(a;)„,,-^(x)„,,|| < (4.8) 

for aWa <^ <1 /a. 

The choice d = a yields (forward) invariance of the cone and the choice d = 1/a 
yields (backward) invariance of the cone C^'^. Note that the cone invariances are in fact 
strict as claimed in the lemma. The above estimates are strict inequalities for x ^ x. 

Now consider the cone in transverse direction, that is '^{x) G C^'^, which amounts 
to \^ {x) c — {x) c\ < cr||\E'(x)tj^s — \E'(a:)u^s||. We have already established invariance. Thus 
|5^c — a;c| < crllxii^s — Xu,s\\ and estimate 04.61) yields 

||^(x),,,-^(x)„,,|| < 

< C-*"-|a;J/'(l + a)||£„,,-x„,,||. 

This is the claimed contraction, K~\ = (^fctum^/Sj^]^ + o") < 1, for 6 small enough. ixi 
Theorem 4.2 The (local) stable set of (0,0,0) under the return map is given by 

"{(x^, Xg, Xc) I Xq Xc(x^, X^)} 

with a Lipschitz continuous function Xc- In particular, Xc(0,Xs) = 0. 

Proof. The idea of the proof is to define a graph transformation on the space of Lip- 
schitz graphs {(x„, Xs, Xc(x„, Xs))}, figure H73l by the inverse return map The cone 
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Figure 4.3: Lipschitz graph Xc = Xc{xs,Xu) with fixed fine Xc{xs,0) = 0. The graph 
transformation given by the inverse return map expands Xu,s and contracts Xc- 

invariance of the previous lemma will ensure that the Lipschitz property of the graph is 
preserved. Due to the expansion/contraction conditions of the previous lemma, the graph 
transformation turns out to be a contraction on the space of graphs. The fixed point of 
this contraction then yields the claim. 

To make this idea precise, consider the Banach space of Lipschitz functions 

X = { C - [0, 5f ^ [-5,5], Xu,s = {Xu,Xs) ^ Xc = C{Xu,s) 

such that C{xs, 0) = 0, Lip(C) < a } 

with sup-norm. The parameters 5,a < 1 correspond to those of lemma 14. 1[ 

Define a map G : X — )■ X as graph transformation, i.e. graph(GC) := \l'~^graph(C). 
More precisely 

GC[{^'\xu,sX{xu,s)))u,s) {^-\x^^sX{xu,s)))c^ for X,,, ^ 0, 

G({0,Xs) := 0. for arbitrary Xs, 

The first equation implicitly assumes that {xu,.s,C{xu,s)) has a pre-image under \l/ and 
that it lies in the domain. The second equation just gives the pre-image of the origin 
under \E'. Note the restriction to non-negative Xu,Xs consistent with the structure of 
the Bianchi system that preserves the signs of the spatial curvature variables Nk- In 
fact due to these constraints the boundary Xg = becomes singular after one excursion. 
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= 0,Xc)u = 0. This fact is not explicitly used in the following but justifies the 
choice of the domain. 

We will prove the following claims 

(i) domain of definition: for all ( & X and Xu,s ^ (0, S] x [0, 6] there exists Xu,.s ^ 
[0,(5]2 \ {0}, such that {<iI-\xu,sX{S:u,s)))u,s = ^u,s- 

(ii) well-definedness: for all C ^ -'^ and Xu,s,S;u,s G [0,5]^ \ {0} the following holds. If 

{^-^{xu,s,C{xu,s)))u,s = {'^~^{^u,s,C{^u,s)))u^s ^ [O'*^]^ ^^^^ already Xu,s = Xu,s- 

Conditions (i) and (ii) yield a well defined function G( for every ( & X. 

(iii) Lipschitz property: for all ^ G X the function G( is again Lipschitz continuous with 
Lip(GC) < (T. 

(iv) contraction: The exists a constant < k < 1 such that for all C, C ^ ^ the estimate 
||G'C-G'C||sup</t||C-Cl|sup holds. 

Conditions (i)-(iii) prove that the graph transformation G indeed maps Lipschitz functions 
in X to Lipschitz functions in X. Condition (iv) provides a contraction. If all four 
conditions hold, then by contraction-mapping theorem there is a unique fixed point, i.e. a 
Lipschitz function C* G X with G(* = (*. Its graph is a locally invariant manifold under 
\E'. It is also the stable set of the origin due to the cone conditions of lemma 14.11 This 
yields the claim of the theorem. Therefore it remains to prove (i)-(iv): 

(i) Let C ^ ^ and Xu,s G (0, 6] x [0, 6] be given. The straight line {xu,s} x [S, C 

is contained in the cone G^^^^Qy We use lemma I4?T1 '${xu,s,0) G S by invariance and 
contraction of the cone Cg'*^. Thus, by invariance and Expansion of C*(x„sO)' image 
of the straight line {xu,s} x [~^7 ^] under \1/ contains a curve in C^(2,^^o) connecting the 
extremal plains {xc = ±5}. By intermediate- value theorem this curve must intersect the 
graph of (. 

(ii) Let C e X and Xu,s,Xu,s e [0,5]^ \ {0} be given with (^"Hx^,^, C(a;„,s)))„^^ = 
{'$~^{xu,s,C{xu,s)))^^s G [0,5]^. Note in particular that {'$~^{xu,s,C{xu,s)))^ ^ 0, as the 
singular line is mapped onto the origin. 
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Thus ^-Hxu,sXixu,s)) e Cl-^, 



C\ 



{^.,.,C(^.,.)) invariance C(5n,s)) G 

The Lipschitz-bound on C G X on the other hand imphes {xu,siC{^u,s)) G 

^(a;„,s,C(x„,s))' "thus {Xu,s-,C{^u,s)) = {Xu,s-, C,{Xu,s)) ■ 

(hi) Again, the Lipschitz-bound on ^ G X translates to {xu,sX{.^u,s)) G C'^xu s c,{xu s)) 
for all x,x. Cone invariance, lemma I^TTl immediately yields the Lipschitz bound on GC,. 

(iv) The singular line is fixed by construction, thus we only have to estimate the 
distance of the nonsingular part. Let C, G X and Xu,si Xu,s £ [0, 5]^ \ {0} be given with 

{^-\xu,sX{xu,s)))u,s = (^-'(5^«,.,C(5^.,.)))„^ e [0,5]'. 

Again, this implies '^~^{xu,s,C{^u,s)) ^ ^^-^{xu s ({xu s)) cone invariance we 

have {xu,sX{xu,s)) e G9 



Thus we can estimate 



^('''tt,*) C('^M,s)l — lie Cll sup "l" II a;^^^ S^Mjsll 

< lie- Cllsup + ^'IC>.,s)-CK.)| 

The first inequality uses the Lipschitz bound on C G X whereas the second one uses the 
cone We obtain 

u , s i'^ u , s ) } 



On the other hand, the expansion of C^-^if^^ 



1-a- 



■IIC-CII 



sup- 



,({Xu,s)) 



under \1/ yields 



(GC-GC) [{^~\Xu,s,aXu,s)))u 

{'^~\Xu,s. C{Xu,s))), - {'^~\Xu,s. C{S:u,s)) 



< 



- Xc(i-(t2)IIC CI 



sup' 



Lemma [4.11 provides constants Kc, a with Kc{l — a ) > 1. Therefore the last estimates 
yields the claimed contraction, k, = l/(A'c(l — o"')), and finishes the proof. [xi 



The stable set given by the theorem determines the initial conditions of trajectories 
that possess the same a-limit dynamics as the 3-periodic heteroclinic cycle: there exists a 
codimension-one set of trajectories that converge to the heteroclinic cycle as they approach 
the big-bang singularity in backwards time. 
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5 Discussion 

We have shown that the set of trajectories within the vacuum Bianchi system, that start 
close to the period-3 heterochnic cycle and follow it in the a-limit t — )■ — oo, form a 
codimension-one Lipschitz- manifold. 

This holds true in Bianchi-IX as well as Bianchi- VIII domains. In fact, the period-3 
heterochnic cycle is a cycle of three pairs of heteroclinic orbits and the domain depends 
on the choice of representatives. 

The same technique is applicable to non-vacuum solutions provided the additional 
eigenvalue fiQ of the Kasner equilibrium is stronger than the unstable eigenvalue, see 
(12. 7p in reversed time direction. In fact the analysis of the local map in section [3] remains 
valid for arbitrary stable dimension and order of stable eigenvalues, provided all stable 
eigenvalues are stronger than the unstable one, just the dimension of the component Xs^ss 
changes. (Note also that the linearization (12. 7p is diagonalizable even at points of the 
Kasner circle at which some of the eigenvalues coincide.) The crucial assumption is a 
positive gap, l/i^l — > 0, between the weakest stable and the unstable eigenvalue of 
the diagonalizable linearization at the Kasner circle. The result on the stable manifold 
of the 3-periodic heteroclinic chain therefore extends to the non-vacuum Bianchi model 
provided 3(2 - 7) > |(1 - ^5), i.e. 

5- 

see (12. 7p and note the reversed time. 

The result also extends to arbitrary periodic heteroclinic chains, as they keep a 
uniform distance from the Taub points Ti, T2, T3. To this end it is important to note that 
all our estimates are uniform on compact pieces of the Kasner circle which do not contain 
the Taub points. Similarly, the Kasner map $ is uniformly expanding on such compacta. 
Therefore, the return map associated to an arbitrary periodic chain is composed of finitely 
many local and global maps as studied in this paper. This yields uniform estimates for 
each of the local and global maps. Hence, the graph transformation of theorem 14.21 can be 
applied to the concatenation of these local return and global excursion maps. As above we 
get a codimension one invariant Lipschitz manifold of points that approach the periodic 
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chain of heteroclinic orbits. Note that these chains may even include heterochnic orbits, 
that approach the Kasner circle in the direction corresponding to the weak eigenvalue. 
The results on the local passage remain valid for in-sections {xg = 1}. We used notation 
consistent with the non-principal direction of the periodic-3 chain, however the order of 
stable eigenvalues was not used in any of the proofs. 

The unstable eigenvalue along the Kasner circle /C in reversed time direction is 
bounded by < 2, see (12. 7p in reversed time direction. The extension of the result on 
arbitrary periodic chains to the non-vacuum Bianchi model is therefore valid as long as 
the eigenvalue /i^ to the eigenvector transverse to the vacuum boundary is stronger than 
2, i.e. 

4 

^<3' 

including all matter models between dust (7 = 1, included) and radiation (7 = 4/3, 
excluded) . 

Non-periodic heteroclinic sequences can be treated in the same way, as long as they 
keep a uniform distance from the Taub points Ti,T2,T3. Thus the limit of the analysis 
presented here is the following result: 

Remark 5.1 Let any heteroclinic chain of Bianchi-II vacuum solutions, i.e. a sequence 
S = ($"(q'))n>o of the Kasner map h2. 5p . he given that keeps a uniform distance from the 
Taub points Ti, T2, T3. 

Then there exists a local codimension- 1 Lipschitz manifold W^IS) of initial conditions 
with backward trajectories in the Bianchi system 112. 1\} following the given chain S . 

More precisely, take such an initial condition, say Xq G W^{S), and its backward 
trajectory {xq ■ t, t < 0}. Then there exists a sequence (tn)n.>0; ^n. < 0, t„ \ —00, such 
that the distance of Xq ■ tn, and $'^(g) goes to zero as n — j- 00. Additionally the pieces 
{xq ■ t, tn+i ^ t < tn} of the trajectory converge to the respective heteroclinic orbits 
connecting $""''^(g) and $"(9), e.g. with respect to Hausdorff distance. 

Note that each transition in the sequence S corresponds to a pair of heteroclinic 
orbits contained in one of the three Bianchi class II ellipsoids. Fixing the sign of the 
nonzero variable Nk in each of the three ellipsoids associates a unique heteroclinic chain 
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to the Kasner sequence. Three equal signs yield a manifold in Bianchi class IX, the 
other choices correspond to Bianchi class VIII. The result holds in all cases and gives 8 
local Lipschitz manifold with boundaries in Bianchi classes VIq, VIIq. These boundaries 
correspond to solutions not following the chain but converging to a Kasner equilibrium 
after one transition. The union of the 8 pieces yields 2 local Lipschitz manifolds without 
boundaries but crossing the subspaces of Bianchi classes VIq, VIIq. The manifolds are 
attached to the 2 heteroclinic orbits given by the first transition of the sequence S. 

Independently, in |BeglO| similar results under additional non-resonance conditions 
of the eigenvalues of the equilibria along the Kasner sequence are obtained by linearization 
techniques. Remarkably, these non-resonance conditions exclude in particular all periodic 
sequences. 

The major question that still remains open in our treatment is the dynamics corre- 
sponding to heteroclinic sequences that come arbitrarily close to the Taub points Ti, T2, T3. 
In fact, generic sequences fall into this last category. Here, the very technical work of Re- 
iterer and Trubowitz |RT10j provides positive results at least in the vacuum case. 

Nevertheless, the question whether the Kasner map already determines the complete 
dynamics of nearby solutions is still open. The answer will probably require a more 
detailed study of the delicate behavior close to the Taub points. 
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